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By explicitly constructing the Hilbert space, Higuchi showed that there is a lower bound on 
the mass of a minimally-coupled free spin-2 field on a curved background [I]. Using the vacuum 
persistence amplitude, we show that this bound is modified by taking into account additional terms 
not prohibited by symmetry in the case of a maximally symmetric spacetime. This result can 
further be generalized to the maximally symmetric space case, such as the FRW universe, and its 
corresponding bound of the deformation parameter is discussed. 


I. INTRODUCTION AND OUTLINE 

The cosmological constant problem is one of the 
most outstanding challenges in theoretical physics. The 
observed value of the vacuum energy density from the 
redshifts of distant galaxies and supernovae is bounded 
from above by 10“‘^®GeV^ (by setting h = c = 1) 
mis, while the change in the vacuum energy density, 
which is expected to occur during electroweak sym¬ 
metry breaking, is ~ 10®GeV^ ; raising the scale via 
renormalization group flow worsens the discrepancy 
such that the predicted change in the vacuum energy 
density becomes ^ lO^^GeV^ at the Planck scale [4]. A 
striking amount of fine-tuning is required in order for 
these scales to be matched with observations. Thus, 
this challenge calls for any kind of a self-consistent UV 
completion of quantum gravity or some modifications to 
either gravity or quantum field theory. In this paper, we 
take modification of gravity (Einstein’s general relavity) 
at linear level (namely the Fierz-Pauli theory), which 
attempts to resolve this issue. A deformation parameter 
(i.e., mass term) to general relativity changes the 
gravitational potential to the Yukawa-type, such that 
its strength decreases exponentially at large distances 
r ~ 1/m, where m is the deformation parameter, and 
thereby suggests a possible explanation to the tiny 
observed value of the cosmological constant in the 
universe [2211^ . 

Since we are interested in phyiscs in the far infrared 
regime and vertices in general relativity (and modifica¬ 
tions thereof), scale as (fc^/Mp^)” [23], it is sufficient at 
first to consider the linearized theory and absorb the 
effects of the interaction in a renormalization of the pa¬ 
rameters. This approximation is implicitly encoded from 
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the long-ranged estimate from the previous paragraph, 
which can be translated into A ~ 1/r ~ m, where the 
deformation parameter m (the graviton mass) is a tiny 
value from the observational constraint comparing to the 
Planck scale jSHSj. A problem arises when one takes the 
massless limit of the linearized theory, called van Dam- 
Veltman-Zakarov (vDVZ) discontinuity [iSl ITO] . which 
causes the discrepancy in certain experimental con¬ 
straints (deflection angle of a light beginning from stars, 
for expample). Since general relativity is a non-linear 
theory, including the higher order terms can smoothly 
recover the massless limit of the spin-2 theory by the 
Vainshtein mechanism [2()j : the massive spin-2 theory is 
only applicable outside the Vainshtein radius, which is 
defined as , where M is the mass 

of static source and Gn is the Newton’s contant. That 
is, the Vainshtein radius sets up a cutoff of the theory, 
and when the massless limit of the graviton is taken, the 
Vainshtein radius becomes infinite, thereby the regime of 
massless spin-2 theory is recovered. Furthermore, there 
are some cases in which the massive spin-2 theory does 
not suffer from the vDVZ discontinuity, such as on the 
negatively curved maximally symmetric spacetime |28] . 

When the degravitation effect on a generic curved 
background is considered, such as the theory of massive 
gravity on a curved background at the far infrared 
regime, one has to pay enough attention to its con¬ 
struction, especially the “covariantization” procedure 
for the Pauli-Fierz theory m- The only safe way of 
covariantizing the Fierz-Pauli theory is to include all 
possible terms that are not prohibited by symmetries and 
keep couplings to terms vanishing in a flat (Minkowski) 
background as free parameters. This can only be 
avoided if the full model including its interactions with 
background geometry is known, which is not the case 
here. Thus, one should consider not only the Fierz-Pauli 
mass term [10], but also terms coupling the field to the 
background curvature. These terms modify theoretical 
bounds on the allowed values of the mass. A possible 
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way of taking this into account is to compute the vacuum 
persistence amplitude (VPA), which is a diagnostic tool 
for testing unitarity of the theory. If the quantum theory 
at the leading order violates the unitarity condition, 
it is unlikely that the full theory would be able to 
cure the violation as far as a legitimate perturbation 
scheme is concerned. We pay attention to de Sitter and 
Friedmann-Robertson-Walker (FRW) spacetimes as the 
given classical backgrounds in this discussion, for those 
are good zeroth order approximations to the universe we 
observe today. 

The program of testing unitarity proceeds as follows. 
We construct the Lagrangian of the linearized massive 
spin-2 field on a classical background in section II. Rie- 
mann normal coordinates are used to obtain a consistent 
approximation scheme. Results are, within this scheme, 
computed up to second adiabatic order. The general ex¬ 
pression of the VPA is given in section III, and by us¬ 
ing the Green’s function derived in the previous section 
along with inserting a localized static external source, 
we evaluate the unitarity of the model. Imposing uni¬ 
tarity gives rise to a new feature in the lower bound of 
the deformation parameter, and its implications are dis¬ 
cussed in section IV. Further analyses and implications 
to another maximally symmetric spacetime with nega¬ 
tive curvature (anti-de Sitter) and maximally symmetric 
space (FRW universe) are discussed in the same section. 
Lastly, more detailed calculations of the Green’s function 
are presented in the appendix section. 


II. LINEARIZED MASSIVE GRAVITY IN A 
CLASSICAL BACKGROUND 


where /i = /i)) = . The bar notation on top of the 

Ricci tensor and scalar implies that these quantities are 
evaluated in terms of the background metric g. Note that 
the last term which carries the coupling constant ^ indi¬ 
cates the coupling between the fluctuation field and the 
background. Since ^ is a free parameter not affected by 
any symmetries, we expect it to be of order one. There 
are several additional ways to couple the Ricci (or Rie- 
mann) tensor to the fluctuation field, but all of them are 
reduced down to this expression for a maximally sym¬ 
metric background. Sm is the mass contribution to the 
action whose expression is given by m 



= -\ J - h^)]. (3) 

In principle, one has to include this deformation term 
with arbitrary relative coefficients for and 

but the only legitimate expression which ensures 
stability of the theory (i.e., absence of ghost) is the given 
expression. Inclusion of the deformation term, along 
with the non-minimal coupling term, seemingly breaks 
the gauge symmetry of the theory (general covarince), 
but this seemingly broken symmetry can be restored by 
the Stiickelberg formalism just as shown in the massive 
spin-1 case of the electromagnetism mm That is, 
there exists a new set of symmetries in the linearized 
massive spin-2 theory by introducing “Stiickelberg 
vector and scalar fields”. The rest of the terms in h with 
cubic orders or higher are included in and our 

interest is only on the first two terms on the right hand 
side of the eq. Q. 


A. Lagrangian of the Linearized Massive Gravity 


The dynamics of the massive spin-2 theory on a curved 
background is obtained from the Fierz-Pauli action [10] 
when we denote the metric of the background by g{x) 
and fluctuation around it with h^^(x). The cosmological 
constant A is included as an explicit source field, and 
the full action is schematically given, after canonically 
normalizing the fluctuation field EZ], by 

S=So + Sm + (1) 


in which Sq contains all quadratic terms of the massless 
spin-2 field 


^0 = 


J d‘^x^/^ 

+l(v-k- 


2V“/iS)Vp/i + 

J - 2A) - - h^) 



Historically, the action of a massive spin-2 field was 
obtained from Einstein’s general relativity. The idea was 
to simply take the Einstein-Hilbert action and expand 
it to leading order in small fluctuations around a given 
background (which in turn solves the full non-linear 
problem). More explicity, the action of the massless 
spin-2 theory on a classical background is derived 
by the perturbative expansion of general relativity, 
g)j.i^ix) = g^i,{x) + h^^(x) along with taking canonical 
normalization of the field. Including the terms up to 
second order in the metric fluctuation yields the 
minimally-coupled (^ = 0) linearized massless spin-2 
theory [5T]. Then, the Fierz-Pauli mass term was added 
by hand to the leading order in fluctuations. From the 
context of effective field theories, this procedure leaves 
a lot to be desired, as this approach does not generate 
the term proportional to ^ in eq. ([^. While the result 
is a consistent theory (at least to leading order) in 
its own, it is not the most general action compatible 
with all symmetries given. Instead, if Einstein’s general 
relativity is to be modified, all allowed terms should be 
taken into account. 

Hence, the full expression of the linearized action for 
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the massive spiii-2 theory at the weak-field limit, with an 
external source added, is summarized as 

S[h;T] = J d'^x^[Co + Cm + h^.T>^% (4) 

and the equations of motion are derived by taking a vari¬ 
ation of the action with respect to the fluctuation field 

■\/ g ~ {^m—^ {hap 9aph) 4” ^ap 

= 0 (5) 

in which /Cm=o is a sum of the kinetic operator and cur¬ 
vature contribution to the massless part. Since the mass¬ 
less spin-2 theoy has two transversal modes, ICm=o is a 
transverse projection operator, and the Bianchi identiy 
(by taking a convariant derivative on the equations of 
motion) imposes four constraints V“/iQ,p = Vp/i with Vp 
being a covariant derivative in the curved background. A 
commutation relation in the covariant derivatives has to 
be found in order to use these constraints, which is 


group, which is orthogonal to the spacelike hypersur¬ 
faces. That is, when the mode expansion is taken on the 
solution to the equations of motion for a quantum field 
in flat spacetime, a basis is chosen such that the solution 
is unique up to a phase because of the timelike Killing 
vector field. As the symmetry of the background is 
changed from Poincare invariance to general covariance, 
however, we do not have such a preferred Killing vector, 
for the manifold is coordinate independent; there exists 
several different Killing equations m ITB] so that the 
notion of “preferred” basis is no longer available on 
a curved background. For instance, when asymptotic 
“in” (for t —oo) and “out” (for t —>■ -l-oo) states are 
defined in the expanding universe, the basis of mode 
expansions for the asymptotic states are non-trivially 
related by the Bogoliubov transformation HU. The vac¬ 
uum expectation value of the particle number operator 
is not necessarily equal to zero if the “in” state number 
operator is evaluated with respect to the “out” vacuum 
and vice versa, thereby a unique particle interpretation 
is lacking [III HU HU. 


[Vct, Va]/lp = TZxahp - Tlxpcah^"’ (6) 

V^Va/lp = VaV^hp - -(5aph - 4hap) (7) 

where the second line is given by specifying de Sitter 
spacetime as the background metric. The equations of 
motion with five independent components are derived by 
using this commutation relation and the constraints given 
above, and taking the trace out (for h is a non-dynamical 
degree of freedom) 


2A 

(—□“ -|- 9^)hap H —^hap 


— Tap + 


1 


3^2 _ 2A L 


A 


^pS/aT - 9 QapT + —gapT ( 8 ) 


in which fj? = + 4^A and = V'^Vp. Let us denote 

this relation by {Oh)ap = (nT)o,p. It is schematically 
solved hy h = G * T {* is the convolution) where G is 
the Green’s function. Since it is difficult to construct 
a Green’s function globally, we compute it only in its 
adiabatic expansion. To do so, first the correct vacuum 
(another difficult issue on general backgrounds) has to be 
identified. This can be done consistently in an adiabatic 
expansion. 


B. Adiabatic Vacuum and Green’s functions 

The concept of a “unique” vacuum state on a generic 
classical background is ambiguous due to general covari¬ 
ance mi uni HI]. Quantum field theory in flat spacetime 
is constructed upon the “unanimous” notion of vacuum: 
A physically unambiguous vacuum state with respect to 
one inertial reference frame is also a vacuum in the other 
inertial reference frames because there exists a preferred 
global timelike Killing vector field dt in the Poincare 


However, an approximate vacuum state can still be de¬ 
fined by imposing certain constraints on the background. 
It is well-known for a slowly varying (or expanding) 
FRW universe that one can solve the equations motion 
by WKB approximation mi HU- In the case of a 
massive scalar field theory in the FRW universe, the 
slowly varying background implies that the period of the 
harmonic oscillator is much smaller than an adiabatic 
parameter, which we call T, and this parameter can be 
chosen as the Hubble time H~^. In other words, we 
take the WKB approximation to the solution based on 
an assumption that the wavelength of the scalar particle 
is much smaller than the curvature scale of the universe, 
Xk ^ and thereby the asymptotic “in” and “out” 

states are defined just as in the flat spacetime. The 
approximation scheme in this paper is up to second 
adiabatic order, and its corresponding vacuum state is 
the second adiabatic one |0(2)) [TT] . 

This approximation scheme becomes more explicit by 
a particular choice of the coordinate system of the back¬ 
ground. We employ Riemann normal coordinate (or a 
local inertial frame), centered at x'. The metric is then 
given by m 

+^T^paup-a{x')y‘^y'’y'" + ■ ■ ■ (9) 

where y = x — x' and the spacetime indices are raised 
and lowered by the flat metric r]p^. We can rewrite the 
expansion scheme shown above as y H~^ ^ A“^/2 (^or 
Aj/2 1) if the background is chosen to be the maxi¬ 
mally symmetric spacetime. In other words, we consider 
a small region around an arbitrary origin x'. This re¬ 
gion is much smaller than the surface of Hubble sphere 
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with radius th ~ and thus is well approximated 

by small corrections to Minkowski. The vacuum struc¬ 
ture, in particular, is very similar to that of fluctuations 
around Minkowski spacetime. Since the curvature scale 
A is a fixed value in the maximally symmetric spacetime, 
the approximation up to second adiabatic order means 
that only the terms up to the quadratic order in y are 
to be included in the action. The de Sitter metric and 
its inverse metric in the normal neighborhood are given 
below 

_ , , A 2 A In 

ga/3{y) = Vay - gVayy + g2/a2//3 + 0{y ) 

rHy) = f 2/"2/^ + 0{y^) (10) 

from which the Christoffel symbol and Riemann tensor 
are computed accordingly. 


With this power counting scheme, we can expand the 
left hand side of the equation of motion up to the second 
adiabatic order as follows 

2A 

(Oh)ap = + y? -\——)hap 

2A 2A 

~ n T M ')^ap T ^X^ap ~^y{(x^ ^p)X 

4A 7 SA^ 2, 2A2 

C'(Q,/lp)A H hap ^^yayph 

4A2 A 7 8A2 ^ 

+ -^2/ 2/(a«p)A - -^Vapy y^hxj 


-^y'^Ohap + ^y^y‘'dpd^h. 


Clip 


( 11 ) 


in which the convention for symmetrization is without 
factor of 1/2. Likewise, the right hand side of the equa¬ 
tion is expanded in the normal neighborhood up to sec¬ 
ond order 


{UT)apiy,x) 


= T, 


1 


dtp 


- 2A 


{( 


h^p^a y Vap T 2 Votp 


T 


rA . A , 2A N 

+ yhpya + ^h^yp ^VapV 


+ 


13”'’ 

/ y'^A _ A^\ 

V^9 


9 

A2s 


dxT 


VapV' 


_ /y?A _ A^x 


V 9 

= (HTo) ap + (HTi) ap + {UT 2 ) ap 


27/ 


2/a2/p 


T 


( 12 ) 


where the last line is defined in such a way that the sub¬ 
script denotes the number of y’s (the adiabatic order). 


Thus, the solution to the equations of motion, evaluated 
up to second adiabatic order, is schematically written as 


hap{x') = J d'^Xx/^{Gyap{x',x){nT)raix) 



1 

+ J2{Gyyp{iiT,y, + (GorypinT^ujis) 

i=0 

and the G{x',x) is the adiabatic Green’s function (sub¬ 
script denoting the order) otbained by solving eq.([TT]). 
Let us perform a Wick rotation x^ ix'^, for it is more 
convinent to compute the VPA in Euclidean space. The 
Green’s function at zeroth adiabatic order is 


/ 74 7 

^ ^ __ Jkix-x') 

(27r)4/2 k^+y^ 

where 5 -^ 1(5^5/- ,5-j;) (14) 


which is the Feynman propagator Ap = l/{k^ + y?) car¬ 
rying indices after imposing constraints {AI°‘hap = ^ph 
as well as taking the trace out) on the equations of mo¬ 
tion in flat spacetime; note that this expression is valid 
only if the effective mass squared is positive, > 0. The 
first adiabatic order Green’s functions are derived using 
the eq. (14) and its defining equations are given as 


[{-D + yyd^JJ]{G^Yx^{y) 

2A 2A 

= -^5i5ydp{G,yyyy) + —y^^Jl^dyGoYxyy) 

-^y"<^(o5p)(G'o)i;(y) (15) 


whose solutions are given by 

{Gyi^k) = l|^(A^^)(Go);^(fc) - ^kyAp)yGyz{k) 

-^.KG^)[Apf - ^i^(-)(A^)3, (16) 

Kjyy'^ = 5;k'^ka + s;k^ka. (i7) 


These are the first order curvature correction terms 
(which is linear in A) to the Green’s functions with re¬ 
spect to the flat spacetime solution. The adiabatic sec¬ 
ond order Green’s functions in momentum space, whose 
defining equations are deduced from the eq. ( [Tl| ), can be 
derived in a similar fashion 


{G2)Zik) = ^-^{ApfiGYZik) + ‘^^Z^VcpiApf - ^ 


81 
108A2 


kyApy{Go)7p{k) 


48A2 

81 


kHApfiGYZik) 


+ KGy)){Apy - ^k^k'^yapiApf - ^y^^kapiAp)* 


81 


81 

+^(iL(7) + KGyykZApf + ^rrfc„fcp(Af )5. 


81 


81 


81 


(18) 
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The detailed calculations can be found in the appendix. 


Therefore, the sum of eq.( 141, (16) and (181 is the inverse 
operator to the left hand side of the equations of motion, 
and it remains for us to put these pieces together into 
the VPA. 


III. VACUUM PERSISTENCE 


Unitarity of the linearized massive spin-2 theory can 
quantitatively be evaluated by using its partition func¬ 
tion whose definition is given |12] 

(0|0)y = (19) 

where h^i, is the massive spin-2 field and its action 
S[h;T] is given by eq.Q, in which T is taken to be a 
classical source; the partition function is to be bounded 
from above by unity if the model is healthy. Since it is 
easier to evalute this quantity in Euclidean spacetime, 
we will take the Wick rotation and see if the partition 
function in Euclidean spacetime is bounded in the same 
way. One might argue that we do not check unitarity as 
the computation will be done in Euclidean spacetime and 
we do not perform the analytic continuation back to real 
times. However, we can still use the result as an indica¬ 
tion of the health of the model. The reasoning is simple: 
We compute a quantity at finite euclidean time. If the 
Hamiltonian is bounded from below (a sign of health of 
the model), the difference between our result and a field 
theory at finite temperature should be small and our 
result should be smaller, but positive. Now, if the ground 
state is stable, at finite temperature the free energy 
F = —ksTlog^Z) (where ks is Boltzmann’s constant, 
T is the temperature, and Z is the partition function) 
should be positive. This implies that Z lies between 
0 and 1. So, our result should be between 0 and 1 as well. 

After taking the Wick rotation, the Jacobian factor 
i /—5 of the integration measure for the action is replaced 
by -i/gE (where E in the subscript denotes the Euclidiean 
signature), and the exponential part picks up another i 
factor along with the action replaced by the Euclidean 
action 


(0|0)r = y’(I?V)e-^-['‘’^l. 


( 20 ) 


This action can schematically be written as 


Se= d'^Xy/gE 




( 21 ) 


where is pairwise symmetric (/i, v and p, a) and 

contains all kinetic, mass and curvature operators. By 
changing the integration variable, 

and after imposing the normalization condition such that 
the amplitude of the linear theory without any external 


source is nomalized to unity, the partition function is 
given below 


(0 I 0)^ = exp{-i J d4xV^[T^-(/C-ir)p,] } 


) exp I 


- / d‘^x^/g]fyx'"''^nu,paX'' 


(0 I 0)° = exp{-i J d4xv^[T^-(/C-ir)p,] } (22) 


in which the second equality is the expression of the 
vacuum persistence amplitude (VPA); the superscript 
0 denotes the normalized partition function. Clearly, 
if a correct normalization is used, the overlap between 
any two states in quantum mechanics of a closed 
system should be bounded by one; otherwise, Born’s 
rule is violated and the overlap cannot be interpreted 
in a probabilistic way. This holds irrespective of any 
manipulations performed with or within the system, so 
the expression in the above constrains the values of the 
exponent, provided the normalization is correct. This 
constraint is equivalent to imposing positivity of the 
Green’s function or absence of tachyonic states in the 
spectrum. 


This general form of VPA becomes more concrete as 
the equations of motion are solved using the power count¬ 
ing scheme Ay^ ^ 1. Inverse operator ]C~^ acting upon 
the external source in the integral is derived by solving 
the equation of motion, and the VPA is replaced by the 
following (the superscript notation of the VPA is dropped 
from now on for the context is clear) 


(0(2) I 0(2))y « exp|-i / d'^x'Jg^d^Xy/gE 

X [r'^"(x')G„(x',x)(nr)^"(a:)] }(23) 


where the integration is limited due to the approximation 
scheme discussed in the previous section and G(x',x) is 
the adiabatic Green’s function up to second order. Let 
us convert this expression to the Fourier space for it is 
easier to calculate the integral. The stress-energy tensor 
at X within the normal neighborhood of x' is 


Tp.y{x) = T^^{y,x) 


d^k 


(2^)4/2^ 


fp^,{k) 


Jk-x' ik-y 


(24) 


where T^^(k) is a Fourier amplitude of T^y{x). The VPA 
in Fourier space representation is, therefore, given order 
by order as follows 


(0|0)r = expj-;^ [ 

^ J\h 


|fe|»Al/2 


d'^k 


0 -|- Gi J- G2]Qp(nTo)T 


+f“^[Go + Gi];^(nfi),<,(A:) 

+r“^(Go);^(nf2) 


( 25 ) 

















6 


where the Jacobian factors, -y/gs and are canceled 

by the inverse Fourier transform of the sources at x and 
x'. Since we are probing the neighborhood region around 
the origin x' whose metric is Minkowskian, the continuity 
equation for the source is reduced to the flat one 


point particle with its mass M localized on a spatial co¬ 
ordinate, (we treat the source as a classical object while 
the spin-2 field is treated quantum mechanically, which 
suffices at the linear level) 


= 0 ^ da^T'^P = 0 


(26) 


and thus kaT°‘P{k) = 0 can be used to simplify the in¬ 
tegrand of eq. (25). The final form of the amplitude is 
given after plugging the eq.(14), (16), (18) and (12) into 
this expression and by specifying the static source as a 


T“^(x) oc (27) 


which apparently obeys the continuity equation so that 
the VPA is given by (detailed calculations are found in 
the appendix) 


(0|0).r =exp{-;^ [ 


-b 


1 

2 J|fc|>Ai/2 
32A2 .4 


a K 

(^ 


+ 


81(3^2 _ 2A) 

20A2 A2 


k^ApY 


r V- 

A /104A 


9(3/r2 - 2A) V 9 


+ 


9(3/r2 - 2A) ^ ^9(3Ai2 - 2Ayr*^^ 

— 16/r^^ k‘^{Ap)^ 

A2 


3 /' 


--(48^^^ - 60A)1 (Apf + , 4 .4-112/r^ + —^^(Af)^ 

81(3Ai2-2A)^ ^ ;jv 81(3^2 _2A)V ^ 3 J ^ ’ 


+ 


48A2(2^2_5a/3),4^^ 

81(3m^-2A) k I 


expressed in the increasing order of the Feynman propa¬ 
gator. It is apparent that the first and second terms of the 
integrand are dominant in the high momentum regime as 
well as the curvature correction terms, which are higher 
orders in the Feynman propagator, do not dominate in 
the infrared regime, for the radius of the normal neigh¬ 
borhood is much sma ller than the curvature radius; i.e., 
k'^ ^ A. From eq. (27), the Fourier amplitude of the 


external source is proportional to a delta function whose 
argument is a zeroth component of the four-momentum 
TP''{k) ~ MJq Jq(5(/co), that the convergence of the VPA 
(i.e., the unitarity) boils down to the behaviour of the 
polynomial function in terms of the Feynman propaga¬ 
tor. Combining the first two terms of the expression gives 


(28) 


r 


the leading order term of the polynomial function as 
ISfP - llA 


9(3Ai2 - 2A) 


(Af)^ 


(29) 


and by substituting k^ with x and with M, factoring 
out the coefficient of the leading order tern and taking the 
common denomimator by the fifth power of the Feynman 
propagator simplifies the expression inside the integral of 
eq. (28) as 


18M-11A fix) 
9(3M - 2A) ' (x + M)^ 

in which expression of the /(x) is given by 


(30) 


fix) = ix + M) 


A 


62 A \ 


32A2 


■ 22M- ^^]ix + M)^ . , 

18M-11AV 3 ’ 9(18M-11A) 

A /104A 108A2M - 100A3, 

+ + + 9(18Af-llA) + 

^ (-112M-b 


x^ix + M) 


+ 


9(18M- llA) 




9(18M- llA) 


(31) 


It is apparent that (x + M)~° in eq. (30) is positive 
because x > 0 and we impose M > 0, so it is necessary 
for us to analyze the behaviour of the polynomial 
function /(x) up to a quartic order in x. In order to 
show the positivity of /(x) in its domain, one has to 


first find its extrema (minima, in particular). Extrema 
are evaluated by taking a first order differentiation 
with respect to x and solving for its roots. Since there 
is only one real extremum of the function, it is the 
point where the minimum occurs (detailed analysis is 
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found in the appendix), given that the coefficient of 
the leading order term, eq. ( [2^ , is positive. It turns 
out that the point at which the minium occurs lies 
outside the domain of the function (x > 0) that the 
actual minimum of the function is at the infrared (IR) 
cutoff of the theory. The IR cutoff is at the scale much 
greater than the curvature scale, and one can take, for 
instance, 3A (or much greater) as the lower bound of the 
integral, and it is shown that /(3A) is strictly positive 
(see appendix). This implies that the integrand of the 
eq.(25) is strictly positive in its domain, and the VPA, 
thereafter, converges to a finite value between zero and 
one. 


If we impose that the VPA does not violate unitarity 
to leading order, we thus obtain a healthy theory if and 
only if 


, 2A 

and 

„ IIA 

(32) 

, 2A 
or M < Y 

and 

, IIA 
0<M <^. 

(33) 

Rewritten in terms of the bare mass m and the 
logical constant, these equations are 

cosmo- 

2 2A 

m + 4^A > — 

O 


> (^ - 4^) A 

(34) 

and 




0 < -1- 4^A < 

18 


-4^A < < (]^- 

-4e)A. 


Depending on the values of certain cosmological 
constants would become unstable. However, there is 
a critical value Co = which is called the conformal 
coupling mm, at which all de Sitter and Minkowski 
spacetimes are stable, as then the hrst condition re¬ 
duces to > 0 (which is the well-known bound from 
Minkowski). 

One further comment is in order: Clearly, due to the 
special choice of the source, it is not unequivocal whether 
the result obtained here actually gives the correct bound. 
There might be a set of unstable states which simply 
do not couple to the source considered here; insofar this 
analysis clearly does not have the claim to completeness 
as Higuchi’s considerations [1]. However, setting the cou¬ 
pling strength to zero, we recover Higuchi’s celebrated 
result, so at least a relevant subset of all “dangerous” 
unstable states couples to source considered here. Our 
bound, in view of this, should be read as a necessary, but 
not a sufficient, condition to obtain stability of spacetime 


under perturbations described by a massive spin-2 field. 
The restriction to a (possibly) small local neighbourhood 
only emphasizes this last point. 

IV. CONCLUSION AND DISCUSSION 


A. Mass Bound 


In his paper [I], Higuchi showed the mass bound of 
the spin-2 theory by the canonical quantization instead 
of the path-integral formulation. A class of equal-time 
commutation relation gives rise to the lower bound 
of the mass squared, > 2/3A, when the negative 

norm state is avoided. However, this lower bound 
physically implies that the quanta of the massive spin-2 
field are able to destabilize de Sitter spacetime, which 
possibly stems from the improper covariantization of 
the linearized theory; a particular choice of the coupling 
between the field and background spacetime, ^ = 0. 
There is no justification of chooing the minimal coupling 
when a general covariant theory in curved background 
is considered. One can recover eq. (341 when the 


non-minimal coupling term is included in the action, 
and perform the equal-time commutation relation just 
as shown in the Higuchi’s paper [1]. The corrected mass 
bound implies the absence of the forbidden range of the 
mass when quantizing the theory. A notion of the global 
vacuum is defined on a spacelike hypersurface in order 
to derive the bound, whereas we used the adiabatic 
Green’s function and adiabatic vacuum up to second 
order of the manifold in order to compute the VPA. So 
it is legitimate if one raise a doubt if the approach in 
this papar is fully valid. Whether the physical set-up 
is global or local, however, as far as any valid vacuum 
to vacuum transition is concerned, eq. (34) is derived 


by imposing the unitarity at the local region, and the 
analysis here is consistent at the linear level. Had the 
theory already broken down in a local region, there 
would be no hope for the same theory to recover globally. 


B. Anti-de Sitter Spacetime 


One can ask whether the anti-de Sitter spacetime is 
also a proper solution the massive spin-2 theory, given 
the de Sitter is its solution, for the previous analysis was 
done in the maximally symmetric 4-dimension spacetime 
regardless of the sign of the curvature. The VPA in the 
anti-de Sitter spacetime is altered simply by replacing A 
with — |A| in eq. (25) as following 






(0(2)|0(2))^ 


= exp 


{-IL 




2^2 + 5|A|/3 


|A| 


,.»|A| (27r)4 L 3^2 + 2|A| 
62|A| 


(Af)- 


4|A| 


+ 


9(3Ai2 + 2|A|) 
20|A|2 


+ 


|A|2 


9(3/^2 + 2|A|) 
32|A|2 




81(3^2+ 2|A|) 


kH^FV- 


|A| 


(-^-16M^)fc^(A^)= 


9(3^2+ 2|A|) 


+ - 


81 81(3/i2 + 2A) 

48|A|2(2m2 + 5|A|/3),4 


(48Ai^+60A) {AfY 


|A|2 


81(3Ai2 + 2|A|) 


(-112/i2 


352|A|\ 


A ^4 

"'f) 


jfc^A 


81(3Ai2 + 2|A|) 


k\AFf]\fik)\^} 


(36) 


and the expression of the effective mass is altered by 
^2 = _ 4^|A|, which implies that the is no longer 

treated as a non-negative parameter. However, the Feyn¬ 
man parameter does not suffer from hitting poles due to 
the domain of the integrand, fc2 ^ |A|, so that we can 
perform the same unitarity analysis on this amplitude as 
in the previous section. The mass bound of the spin-2 
field comes from the positivity of the coefficient of the 
leading order term 

18/r^ + ll|A| ^ 2 . 11|A| 

9(3^2+ 2|A|) ^ > 18 

m2>(4e-^)|A| (37) 

and the lower bound to the mass is ruled out by the 
similar reasoning given the previous section. This mass 
bound looks again non-trivially related to the background 
curvature when the coupling strength is set to equal 
to one. We can constrain this value within a range 
— 1 < ^ < 11/72 ~ 1/7, and the mass squared is again 
bounded from below by a negative quantity. Notice that 
the massive spin-2 field at the linear level is weakly cou¬ 
pled to the background geometry and that, if the coupling 
strength is adjusted in the de Sitter spacetime, the same 
theory is to be held in the anti-de Sitter spacetime with 
the similar coupling constant, given that both are proper 
solutions to the massive gravity. This compatibility be¬ 
tween the two maximally symmetric spacetimes suggests 
the restriction on the value of the coupling between 1 /6 
and 1/7 in such a way that the constant becomes slightly 
greater than (or equal to) 1/6 when the background has 
a positive curvature and it is slightly smaller than 1/7 
if the background has a negative curvature. Imagine 
a free massive graviton is propagating in the de Sitter 
background, and all of a sudden, the background is lo¬ 
cally changed into the anti-de Sitter spacetime due to the 
presence of a source, such as a black hole m- The same 
graviton is also to be freely propagating in this back¬ 
ground regardless of its curvature, and these two bounds 
of the coupling constant ensure the safety of the graviton 
in the maximally symmetric spacetimes. Had the mass 
had a certain forbidden region, as is claimed in [1], the 
quantum description of the graviton with a fixed value of 
the mass is to be broken down by the unitarity violation. 


when it enters another region of the constant curvature 
whose magnitude is greater than the graviton mass. This 
violation implies that either the background is unstable, 
the degrees of freedom are altered by the change of the 
background, inclusion of the higher order terms might 
restore the unitarity, or the unitarity can be restored in 
some other ways. In any case, the effective description 
is to be modified in some way or the other when the 
unitarity violation occurs. 

C. FRW Universe 

The unitarity analysis given above is slightly modified 
when the background is generalized from maximal sym¬ 
metry in spacetime merely to space, such as a homo¬ 
geneous and isotropic universe (Friedmann-Robertson- 
Walker or FRW universe in short). Recall the relation 
between the normal coordinates (or the wavelength of 
graviton) y and the Hubble time H, for a time parame¬ 
ter in FRW universe is y'^ <C H~^. The Hubble param¬ 
eter is a constant value in de Sitter and anti-de Sitter 
spacetimes, but it becomes a time-dependent parameter 
in the comoving frame once the spacetime symmetry is 
abandoned; i.e., H‘^{t)y'^ ^ 1 is the adiabatic parameter. 
When the universe is assumed to be spatially flat, which 
is what we are observing today pHEil] , and its expansion 
rate is slow {H <C H'^), the adiabatic parameter can also 
be written in terms of the curvature scalar 

7^(^)y2 < 1 (38) 

which clearly differs from de Sitter (or anti-de Sitter) case 
such that the Ricci scalar is now time dependent. Slow 
expansion {H ^ H^) of the universe implies that the 
significant expansion takes place in the time span much 
larger than the Hubble time), so that we can mimic the 
similar analysis in the maximally symmetric spacetimes, 
for the curvature scalar is almost a constant in the normal 
neighborhood region. Therefore, the mass bound in FRW 
universe naturally arises when unitarity is imposed 

m^> (39) 

and this agrees with the literature Dni], when the cou¬ 
pling strength is minimal. This bound seemingly runs 
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due to the time dependence of the Hubble parameter. 
We can treat the Hubble parameter as a constant in this 
context if the local region of the spacetime is restricted 
enough, and conclude that the mass bound is again ab¬ 
sent in the slowly expanding and spatially flat FRW uni¬ 
verse. 


V. APPENDIX 


A. Adiabatic Green’s function 


As taking the Fourier tansform of the eq. (151, the first 
term on the right hand side of the equation is expressed 


as following 


2A 




where 


8 A 

T 


^ap 




- 


4A 

T 


(fc 2 -I- ^ 2)2 


8^ = ^ 
dk0 


(40) 


and Sm is defined in the eq. (14). The second and the 


■^ctp 


third terms are transformed in a similar fashion 


- —o. 


(a 


STcr 

q °A|p) 

" A: 2 -P^ 2 J 

9 Lfc 2_|_^2 


4A 


gTO- 


K 


(ro-) 


9 Lfc 2 -|- _|_ ^ 2^2 


^OLp 


k: 


Otp 


8 A - 

9 ~k'^ + fj? (k^ + 
(t 


(41) 


in which the definition of Kpa ' is given in eq. (16), 


and thus, the first order Green’s functions are derived 
as in the eq. (16). After solving for the first order 


Green’s function, the defining equations of the second 
order Green’s functions are given as 


+ = - — 6^6y^^iG^)l’:;iy) + —y^^6;^d\G,)i;iy) - —y^6jMG^)V,iy) 

042 94 A^ 

-^y^st6;{Go)7,{y) + -^ycypri^^Go)7,{y) - -^y\c.s;^iGo)7,{y) 

8A2 2A a 

+^Vapy^y^Go)7^{y) + —y^a6:6'^^{Go)7,{y) - -y^y-^dpdJl5;{G,)Y,{y). (42) 


Notice that the last two terms on the right hand side of the equations cancel each other because the zeroth order 
Green’s functions are Lorentz invariant; i.e., they are proportional to S^^y^ that both terms are identical with the 
opposite sign when we plug this expression into the equations. The first term on the right hand side of the equation 
in Fourier space in terms of the zeroth order Green’s functions is given by (let us keep 1/81 as a common numerical 
factor in front of each term from now on) 


9 A QfiA^ 144A^ 4SA^ 

d^[kp{G^)i;)m = —(A^)(Go);^(fc) - -^k\AF)HGo)Z{k) + —fc"(Ap)3(Go);^(fc) 


81 
24A2 

81 


81 


81 




24A2 

81 




(43) 


and likewise the second and third terms are transformed as 


9A 4RA^ ‘V) 79A^ 

- = - — {Ap){Go)Z{k) + —k\Apf{G,)Z{k) + —{kGZ+kGz))[AZ^ 


81 ^ ^ 81 
16A2 48A2 

-k k rjapiAp) - 


81 


96A2 . 


(^( 7 ) + KGz))kZApf - —k-k^k^kp{Apf (44) 


pOi 


81 


81 ' ' 81 
4A 96A2 32A2 172A2 

dHUkp){Gi)Zm = —iAZ{Go)Zik) - ^iAZ\Go)Z{k) - 


81 
16A2 


96A2 


81 


Z’^k^kpiApf + —{kZZ^ + KZZ^)kZAF)^ + 


192A2 


rfc'"fc„fcp(AF)'‘.(45) 


81 


81 


81 
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At last, the remaining terms (from the fourth to the seventh) are given by 

8 A^ 

2 A 2 

81 + 


-^{AF){Go)Z{k) + ^k^iAFriGo)Z{k) 


4A2 

~W 


rrcT 


(Af )(G„)S(t) + 


rrcr 

°A7 


fc 2 + /J ,2 


18 

16A2 


81 


81 


(\A\2 

- —r^p^pk^k^Apf 


(46) 

(47) 

(48) 

(49) 


from which one can solve for the second order Green’s function by summing eq. (43), (44), (45), (46), (47), (48), and 
(49) and divide it by the Feynman propagator, whose expression is given in eq. (18). 


B. VPA expansion 


Let us calculate the term by term in the integral from the eq. (25) with the given Green’s function (eq. (14), eq. 


(16) and eq. ([l^). Referring to the eq. (12), the last term from the eq. (25) in the Fourier space is given by 


f‘^^{k){G^)Zm2)ra{-k) = 


- A 2 

27(3^2 _ 2 A) 


r“^[277„p(Af)2 - Srj^peiApf + 8fc„fcp(A;^)3]f (-fc) 




(50) 


in which it is clear that the conservation of the source kaT°'^{k) = 0 is used on the second line. One can perform the 
similar calculation on the second term in the eq. (25) as follows 

f“^(fc)[(Go) + {G,)]ZiIlf,)rA-k) = ,, [^fc^(A^)^ + 


k^iAFf-^^k\AFy]\fik)\^. 


The first term in the eq. 

f‘^Zk){Go)Z,i^To)rA-k) = - 


3ii2-2AL 9 ^ ^ ' 81 ^ 81 

is given order by order such that the zeroth order term is given by 


(51) 


1 


(A^^-f)|rp](A^) 


3/r2 - 2A 


2ii2 _ 5A/3 
3/r2 - 2A 


(A^)|f(fc)|2 


(52) 


and the expression of the localized source (eq. (27) was used for the second equality. Likewise, the first order Green’s 
function contribution to it is 


r“"(fc)(Gi);"(nfo)..(-fc) = 


12A(2/i2 - 5A/3), ^ 2 4A(2ii2_5A/3) 2 .. ^3 




kZAZ^ |T(fc)| 


(53) 


9(3/i2-2A) ' 9(3/r2-2A) 

where the identity KaZ^T°‘^{k) = 0 from the conservation equation is used. The second order contribution is the 
following 


r“"(fc)(G 2 );^(nfo)..(-fc) = { 


20 A 2 A 2 (48/^2 - 60A) 


L 81 


81(3/r2 - 2A) 


(Af)= 


A 2 


81(3/i2 - 2A) 


(^—k^iAp)^ + k^{AF)^^ 


+ 


A 2 




81(3Ai2 - 2A) 
is altogether given by 


+ 


352A 


fc"(AF)" + 


48A2(2^2 _ 5^/3) 

81(3Ai2 - 2A) 


k^Apf 


'(54) 


so that the first term of the eq. 

f“GMrr ir iri-rnri r ^-fV-sA/d 12 A ( 2 m 2 _ 5^/3 4A(2/.2_ioa/9) 2 ., .3 

T (fc)(Go + Gr+G.)„,(nro)..(-fc)-|^^^^-^(A^) + ^^^-^-^(A^) " 9 ( 3^2 _ 2A) ^ 


r20A2 A2(48/i2 - 60A) 

81 81(3^2 - 2A) 

48A2 


3 A2(-112/r2 + 352A/3) 2 ., ^4 


+ 


' 81(3/i2 - 2A) 

48A2(2/x2-5A/3). 


Lfc 2 (A, 


81(3/i2 - 2A) 




(55) 


and the VPA in the momentum space is written in the 


C. Bound Analysis 


eq. (28). 


As the quartic function f{x) is expressed in eq. (31) 
and its positivity on its domain are concerned, its ex- 
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trema are found by taking first order differentiation with 
respect to x and setting it equal to zero (can be checked 
by Mathematica). Two extrema are complex values, and 
there only remains only one real extremum, xq, which 
is the point at which the minimum occurs. Since the 
condition 


18M- llA 
9(3M - 2A) 


(56) 


is assumed from the beginning, two different plots of 
the minima Xq with respect to the two parameters, 
A and M, are to be analyzed: The contour map for 


is positive in its domain when the two conditions of the 
parameters are assumed. 




0 1 2 3 4 5 

A 


FIG. 1. A(= A) versus M when M > 2A/3. Colored regions 
signify the minima xq, with the contour lines progress as -2, 
-4, and so forth; i.e., the evaluated minima are out of the 
domain of the function f(x). 


FIG. 3. /(3A) when M > 2A/3. Colored regions are the val¬ 
ues of /(3A), and the contour lines progress as 25000, 50000, 
and so forth, which means f(x) is strictly positive in its do¬ 
main, given M > 2A/3. 

One can choose an arbitrary value of x as long as it 
is much greater than the curvature constant A (recall 
X ^ A from the adiabatic expansion), so let us choose 
3A as a test. A contour plot for /(3^ in terms of A and 
M for M > 2A/3 is given in FIG. [^such that the first 
contour line corresponds to /(3A) = 25000, the second 
line to 50,000, and so forth; f{x) is strictly positive for 
X = 3A. 


M > 2A/3 is given in FIG. in which the horizontal 
axis A is the value of the curvature constant normalized 
by the corresponding unit and vertical axis is the value 
of M normalized likewise. There is a forbidden region 
between the colored ones, and this is where xg is not 
well-defined from the condition f'{x) = 0. Contour line 
at the bottom corresponds to —2, the next line to —4 
and so forth, from which we conclude Xg < 0. The 
second plot for the other condition, M < IIA/18, has 
a simpler feature where its contour line is at —2 and 
the rest are also negative values (see FIG. [^. These 

6 
5 
4 

S 3 
2 
1 
0 

0 1 2 3 4 5 

A 

FIG. 2. A(= A) versus M when M < llA/18. Golored 
regions are the minima xo, and the contour line corresponds 
to Xo = —2. The rest are also negative values, which means 
the minima is out of domain of the function f(x). 

two contour maps inform us that the minimum of f(x) 
occurs when Xg < 0, which is outside the domain of the 
function, so it is necessary for us to check if the function 


6 
5 
4 

S 3 
2 
1 
0 

0 1 2 3 4 5 

A 

FIG. 4. /(3A) when M < llA/18. As in FIG. 3., colored 
regions correspond to the values of /(3A), and the contour 
lines are given as 25000, 50000, and so on. In other words, 
f(x) is strictly postivie in its domain in this case, as well. 

The second contour plot (FIG. shows a similar 
feature such that its first contour line corresponds to 
/(3A) = 25, 000, the next line is at 50,000, and so on. 
One can go on with greater values of x, and it shows 
the similar behaviour. Therefore, we conclude that the 
function f(x) is strictly positive in its domain, which is 
sufficient to ensure the unitarity of the VPA in the lin¬ 
earized massive gravity. 
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